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We provide a thorough exposition of recent results on the quantum stabiUzation of cosmic 
strings. Stabilization occurs through the coupling to a heavy fermion doublet in a reduced 
version of the standard model. The study combines the vacuum polarization energy of 
fermion zero-point fluctuations and the binding energy of occupied energy levels, which 
are of the same order in a semi-classical expansion. Populating these bound states assigns 
a charge to the string. Strings carrying fermion charge become stable if the Higgs and 
gauge fields are coupled to a fermion that is less than twice as heavy as the top quark. 
The vacuum remains stable in the model, because neutral strings are not energetically 
favored. These findings suggest that extraordinarily large fermion masses or unrealistic 
couplings are not required to bind a cosmic string in the standard model. 
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1. Introduction 

It is well-known that the electroweak standard model and many of its extensions 
have the potential to support string-like configurations that are the particle physics 
analogs of vortices or magnetic flux tubes in condensed matter physics. Such objects 
are usually called cosmic strings to distinguish them from the fundamental variables 
in string theory, and also to indicate that they typically stretch over cosmic length 
scales. They are also called electroweak strings or Z-stringJIHSI because the Z- 
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component of the electroweak gauge boson acquires the structure of the Abehan 
Nielsen-Olesen vortejPl. 

Such strings may have emerged copiously in the early universe at interfaces 
between regions of different vacuum expectation values of the Higgs field(s) in sev- 
eral stages after the Big Bang until electroweak symmetry breaking was reached. 
If strings are absolutely stable they should have survived and we should be able to 
observe them today. In the electroweak standard model topologically stable strings 
are ruled out, but the absence of this particular stabilizing mechanism does not 
imply that electroweak strings are unstable or irrelevant for particle physics. While 
their direct gravitational effects ar e ne gligible, Z-strings can still be relevant for 
cosmology at a sub-dominant levei^^. Their most interesting consequences origi- 
nate, however, from their coupling to the standard model fields. Z-strings provide a 

o 

source for primordial magnetic fields ^ an d they also offer a scenario for baryogenesis 
with a second order phase transitiorl^^^. In contrast, a strong first order transition 
as required by the usual bubble nucleation scenario is unlikely in the electroweak 
standard model^ without non-standard additions such as supersymmetry or higher- 
dimensional operator JI^. 

These interesting effects are only viable if Z-strings are energetically stabilized 
by their coupling to the remaining quantum fields. The most important contri- 
butions are expected to come from (heavy) fermions, since their quantum energy 
dominates in the limit Nc oo, where Nc is the number of QCD colors or of 
any other internal degree of freedom. The Dirac spectrum in typical string back- 
grounds is deformed to contain either an exact or near zero mode, so that fermions 
can substantially lower their energy by binding to the string. This binding effect 
can overcome the classical energy required to form the string background. How- 
ever, the remaining spectrum of modes is also deformed and for consistency its 
contribution (the vacuum polarization energy) must be taken into account as well. 
Heavier fermions are expected to provide more binding since the energy gain per 
fermion c harg e is higher; a similar conclusion can also be obtained from decoupling 
arguments^. 

A number of previous studies have investigated quantum properties of string 
configurations. Naculicfff2] has shown that in the limit of weak coupling, fermion 
fluctuations destabilize the string. The quantum properties of Z-strings have also 
been connected to non-perturbative anomalieJI^. The emergence or absence of ex- 
act neutrino zero modes in a Z-string background and the possible consequences 
for the string topology were investigated in Ref. IT. A first attempt at a full cal- 
culation of the fermionic quantum corrections to the Z-string energy was carried 
out in Ref. ,15^ Those authors were only able to compare the energies of two string 
configurations, rather than comparing a single string configuration to the vacuum 
because of limitations arising from the non-trivial behavior at spatial infinity. We 
will discuss a solution to this obstacle in detail below. The fermionic vacuum polar- 
ization energy of the Abelian Nielsen-Olesen vortex has been estimated in Ref. [T6| 
with regularization limited to the subtraction of the divergences in the heat-kernel 
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expansion. Quantum energies of bosonic fluctuations in string backgrounds were 
calculated in Ref. |17i Finally, the dynamical fields coupled to the string can also 
result in (Abelian or non-Abelian) currents running along the core of the string. 
The time evolution of such structured strings was studied in Ref. 1181 where the 

current was induced by the coupling to an extra scalar fi eld. 

This presentation is based on a number of publication|I^H22! regarding the con- 
tribution of fermion quantum corrections to the vacuum polarization energy of a 
straight and infinitely long cosmic string that appeared over the last year or so. 
Technical details underlying the results given here can be learned from those pub- 
lications, cf. in particular the appendices of Ref. [22l 



2. Model 

For the current investigation the fermion doublet will be considered degenerate so 
that the introduction of a matrix notation for the Higgs field is appropriate. Then 
the string configuration reads 

-icosi^i) sin(ei)e™^y ^ ^ 

for the Higgs field and 

M/ ■/,jGip)^f sin(6) zcos(6)e-"''^\ 

for the gauge boson (in temporal gauge) . The variables p and ip are polar coordinates 
in the plane perpendicular to the string, while the Higgs vacuum expectation value v 
and the gauge coupling constant g are model parameters. The string configuration 
involves profile functions Jh and fc which are the analogs of the Nielsen-Olesen 
vortex profiles with boundary conditions 

P — ^ : fa, Ih — ^ 

P — ^ oo : fa, Ih — ^ 1 • 

The integer n is the winding of the string, for which we will typically take n = 1 in 
numerical calculations. Finally, the (variational) parameter measures the relative 
strength of the Higgs and gauge boson components of the string. 

We will consider a modified version of the electroweak standard model which 
only has a SU(2)l gauge symmetry. This modification is equivalent to a vanishing 
Weinberg angle. Then the (classical) boson part of the Lagrangian reads 

C^os = -^tr (G^'^G^,) + itr (I?'^<f>)^ D^'f - ^tr (a>t$ - v^f , (4) 
with the covariant derivative Z?^ = 9^ — igW^ and the SU{2) field strength tensor 



(5) 
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The boson masses are determined from g and v and the Higgs self-couphng A as 
mw = gv/\/2 and mn = 2t) -s/A for the gauge and Higgs bosons, respectively. The 
interaction with the fermion doublet is described by the Lagrangian 

Ae, = i^{PlP+ PrP) * - / * + * with Pr^l = i (1 ± 75) . (6) 

Upon spontaneous symmetry breaking the Yukawa coupling / induces a fermion 
mass m ~ vf. The similarity with the standard model of particle physics suggests 
the model parameters 

5 = 0.72, 7;=177GeV, mH = 140GeV, / = 0.99 , (7) 

where we have taken the fermion doublet to have the mass of the top quark. In our 
numerical search for a stable string, we will also study other model parameters, cf. 
section [5] 

The classical energy per unit length of the string is solely governed by £bos: 



— 2t: I pdp Iri^ sin^ ^1 







9 \ P J J P 



-f + (8) 



where the dimensionless radial integration variable is related to the physical radius 
by Pphys = p/m and p.H = mii/m. 

The central object of the present investigation is the fermion contribution to 
the energy. It will obtained from the solutions to the Dirac equatior@ in the two- 
dimensional plane perpendicular to the string 



iJvf™ = w„vl,„ with H^~i(J^/'P]dp--(^^.'''^'^]d^ + Hint, (9) 

'fT-p J p \(T ■ If ' 



where the single particle Hamiltonian, H is extracted from Cfcr- The interaction 
part, iJint depends on the chosen gauge. A specific choice of gauge, which accom- 
modates the subtleties at large p, will be discussed in Sec. 13.31 The spectrum of the 
Dirac operator consists of bound state solutions with discrete eigenvalues cj„ = ej 
and continuous scattering solutions whose eigenvalues uj are labeled by momentum 
fc, i.e. Lu = \/W^-\-rn? . Finally, the Dirac Hamiltonian for the string background 
anti-commutes with — 7*^7'^. Thus the spectrum is charge conjugation invariant 
and it suffices to consider the non-negative eigenvalues. 



''We use the standard representation with 70 = diag(l, 1,-1,-1). 
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3. Formalism 

The vacuum polarization energy is the renormahzed sum of the changes of the 
zero-point energies of fermions in the background of a static configuration 



E-var. — 



-IE 



9 2^ 



(10) 



Here Un are the energy eigenvalues in the presence of the string as obtained from the 
Dirac equation and wi^'' are their free counterparts. The overall sign has been 
chosen to describe fermionic vacuum fluctuations; for bosons it must be altered. In 
the second part of equation ([TU]) the changes of the single particle energies are re- 
written as the contribution from distinct bound states (e^ ) and the (renormahzed) 
change of the density of scattering states, A Prcn{k). These states obey the standard 
dispersion relation ujk = + rn? with the mass m of the quantum fluctuations. 
In the above, the factor h has been made explicit to stress that -Evac is a quantum 
effect. As in the previous section natural units [h ~ 1 and c = 1) will be adopted 
in the remainder of this paper. 

The computation of the fermion contribution to the vacuum polarization energy 
of a cosmic string proceeds in three stages. First, spectral methods are employed 
to express this energy in form of scattering data^. The equivalence of terms in 
the Born series and Feynman diagrams makes it possible to impose standard renor- 
malization conditions. In a second step this approach is extended to accommodate 
configurations that are translationally invariant in one or more spatial directions, 
using the so-called interface formalisrrPSl. Third, we have to cope with the fact 
that the cosmic string configuration has a non-trivial structure at spatial infinity. 
Because the fields approach a pure gauge rather than zero, the direct application 
of spectral methods is impossible. Rather a particular subset of gauges has to be 
adopted to define a well-behaved scattering problerrPQl. 



3.1. Spectral Methods 

The basic idea of the spectral approach is to express A /9ren(^) as the momentum 
derivative of the phase shifts of the scattering stateJ^. Then 

e..c^~^Y.',~J2d^I ^^^d^t-^^^^)]-- (11) 

It is implicitly assumed that the system allows for a partial wave decomposition with 
degeneracy factor Di for the partial wave of (generalized) angular momentum £. 
The major concern in eq. (|lip is renormalization. As it stands, in three spatial 
dimensions the vacuum polarization energy, eq. (jlip . is quadratically divergent in 
the ultra-violet. For large momenta the Born series adequately represents the phase 
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shif10. We therefore consider 

^^"^ - E - E ^ i I Mk)]N , (12) 

where the subscript on the phase shift indicates the subtraction of the N leading 
terms of the Born series to the phase shift: [5e]j^ — 5i — 5^p — 5^ — ... — (5^^'. 
Choosing TV sufficiently large renders E'^^^ finite. Technically the Born series is an 
expansion in powers of -ffint- The contribution of any such power to the vacuum 
polarization energy can be associated with a Feynman diagram 



where the loop corresponds to the quantum fluctuation and the external lines repre- 
sent insertions of the background field, i.e. the interaction of the cosmic string with 
the fermions. Therefore, the artificial subtraction in equation (jl2p is the sum over 
all Feynman diagrams with up to N insertions of the background interaction: -Bp^' . 
These diagrams can be computed with standard techniques, in particular dimen- 
sional regularization can be implemented to handle the ultra-violet divergences. 

The most important advantage of shuffling the divergences into Feynman di- 
agrams is, however, that they can be straightforwardly combined with the coun- 
terterm contribution to the energy, i?cT- The latter is found by substituting the 
background configuration into the counterterm Lagrangian in just the same way as 
the classical energy, eq. ([8]) is obtained from Chos ■ In any (multiplicatively) renormal- 
izable theory the counterterm Lagrangian has the same structure as the classical 
one and a suitable choice of its coupling constants cancels all ultra-violet diver- 
gences completely. Then E^ + Eqt is free of any divergences. The finite pieces of 
the counterterm coupling constants will be uniquely determined from appropriate 
conditions describing properties of the particles that are associated with the fields. 
This procedure will be briefiy discussed in section. 14.31 In total, the sum 

i?vac = E^''^ + + EcT (14) 

gives an unambiguous result for the vacuum polarization energy once the renor- 
malization conditions are fixed. A first principle derivation of this result based on 
a quantum field theoretic formulation of the energy momentum momentum tensor 
rather and the analytic properties of the Greens function for scattering boundary 
conditions is presented in Ref. [26l 



) ■ , (13) 



'^For small momenta it does generally not converge towards the exact phase shift, in particular 
when bound states are present. 
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3.2. Interface Formalism 



Equation (jl4p is perfectly suited for a computation of the vacuum polarization 
energy for static configurations that allow for a full decomposition into partial waves. 
However, this is not the case for the cosmic string which is translationally invariant 
along the z-axis in coordinate space. In this scenario the wave-function of the 
quantum fluctuation factorizes into 



^ix,t)^e~'^'e'P'Mp) (15) 

so that we have a scattering problem in the plane perpendicular to the string, and 
the total dispersion relation a; 

replacements ej 



This implies the 



+ ^1 in equation (jl2p so that an intc- 



ip" + e| and LJk ^ y y ^ "^j, 
gration overp (with measure dp /2tt) yields the vacuum polarization energy per unit 
length of the string. In doing so an immediate obstacle arises. The scattering data 
do not depend on the momentum p along the symmetry axis. Hence for any A^, this 
p-integral will not be finite. A careful analysis treats the p-integral in dimensional 
regularizatior 



TT ^ 



(16) 

where d is the analytic dimension of the subspace in which the configuration is 

translationally invariant. The divergence now manifests itself via the sin gula rity of 

1971 

the F-function coefficient as o? 1 . Due to sum rules for scattering data^, which 
represent generalizations of Levinson's theorem, the expression in curly brackets in 
eq. (fT6|) vanishes as d ^ 1. Hence this limit can indeed be takeri^'^. Finally, the 
fermion spectrum in the string background is charge conjugation invariant. Thus 
i?^^^ is twice its contribution from the non-negative part of the spectrum. 



\ f"° dk 

















dk 



Mk)] 



N 



■E 



' In-^ - e2 



(17) 



Here is an arbitrary renormalization scale that has no effect on E'^^^ by exactly 
the same sum rules. The expression, eq. (|17p . for E^^'^ replaces the analog in equa- 
tion p4| . Note that the function multiplying the (Born subtracted) phase shift is 
of higher power in k than its counterpart before integrating over the momentum p 
conjugate to the coordinate of translational invariance. Hence N must be increased 
if directions are added in which the configuration is translationally invariant. This, 
of course, merely reflects the fact that ultra-violet divergences turn more severe in 
higher dimensions. 

We stress that the expressions obtained so far for the vacuum polarization en- 
ergy strongly rely on the analytic properties of the scattering data. Furthermore, 
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in numerical calculations, we mostly consider them for purely imaginary momenta. 
This has (at least) two advantages: (i) the oscillating phase shifts turn into ex- 
ponentially decaying (logarithms of the) Jost functions, and (ii) the momentum 
integral and the sum over angular momenta may be exchanged. While (i) drasti- 
cally improves numerical stability, (ii) significantly simplifies the treatment of the 
logarithmic divergences that emerge at third and fourth order of the Born and 
Feynman expansions. These technical details are discussed at length in Ref. .22 and 
briefly addressed in section 13.41 



3.3. Choice of Gauge 



We have now established a formalism for computing the vacuum polarization en- 
ergy of background fields in the string geometry. However, we still have the prob- 
lem that the string does not induce a well-behaved scattering problem because 
of its non-trivial structure at spatial infinity. Though gauge invariant combina- 
tions of the Higgs and gauge bosons are trivial at spatial infinity, the individ- 
ual terms in the Born and Feynman series are not gauge invariant and there- 
fore ill-defined. This ambiguity appears because the Dirac Hamiltonian that is ob- 
tained by straightforward substitution of the field configuration, eqs. ([I]) and ([2]) 
does not turn into the free Dirac Hamiltonian as p — )■ oo, but instead becomes 
H — ^ W{(p)HhccU {(f). This local gauge transformation acts only on the left-handed 
fermions, U{ip) = P^exp {in ■ t£i) -I- Pr with fi = {cos{nip), — sin(n(/j), 0) . Unfortu- 
nately, the gauge transformation H — U{tf)HU'^{if) does not solve the problem 
for all p e [0, oo]: Although it would generate vanishing interactions at infinity, it 
also induces a l/p^ potential at the center of the string, p Q. This might still 
yield well-defined phase shifts, but the conditions underlying the analyticity of the 
scattering data are certainly violated by this singular behavior. As argued at the 
end of the previous section, analyticity is central for numerical feasibility of our 
approach. As a solution, we can define a radially extended gauge transformation 

U{p, ^) - PiCXp {ih ■ Ti{p)) + Pr . (18) 

it yields the interaction 



This transformation fixes the gauge and in equation 
term 



cos(A) 



]| 

-I 



]| 

-I 



2 dp 



a 



a ■ p 



P 

a-p 



n 

Yp 



—a ■ Lp 
G ■ (p 



a ■ 

-a 



i sin(A) 

/Gsin(A)/G(A)-f (/G-l)sin(0/G(-e) • (19) 



The new gauge function ^(p) is hidden in the difference A(p) = £i — £(p) which 
appears both explicitly and as the argument of the space-dependent weak isospin 
matrix 

^ — sin(.T) — i cos(a;) e*"'^\ 
cos(x) e~*"'^ sin(x) j 



Ig{x) 



(20) 
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Imposing the boundary conditions ^(0) = and ^(oo) = for the new gauge 
function ^(p) defines a well-behaved scattering problem. Otherwise, the specific 
form of ^(p) is irrelevant. This property allows us to verify our numerical results by 
modifying its shape while keeping the boundary conditions fixed. 

All explicit matrices in eq. (jl9p act in spinor space. Together with the bound- 
ary conditions, eq. a well-behaved scattering problem is obtained. With this 
choice of gauge a scattering matrix and, more generally, a Jost function can be 
straightforwardly computed. Moreover, the Born series to these scattering data can 
be constructed simply by iterating Hmt- 

Note that the gauge transformation is single-valued at spatial infinity, 
U{oo,ip) = U{oo^ip + 27r). In this respect it differs from the analogous problem 
of fractional fluxes in QED. In that case a similar choice of gauge is hence not a 
remedy; rather the calculation of the vacuum polarization energy require s th e in- 
troduction of a return flux to arrive at a well-behaved scattering problerrPH. The 
return-flux approach can also be used for the present calculation, but it is much 
more laborious numericalljff^ISni 



3.4. Fake Boson Field 

The idea of utilizing a fake boson field to simplify the treatment of higher order 
divergences was first implemented in Ref. 29, As mentioned above, the continuation 
to imaginary momenta k ^ it and 5f{k) — > where is the logarithm of 

the Jost function, allows the exchange of the momentum integral with the angular 
momentum sum. Then the third and fourth order contribution from the Born series 
produce logarithmic divergences. These divergences are similar to the ones found 
in the second order vacuum polarization energy of a boson field fluctuating about 
a scalar potential. Matching its strength appropriately allows to replace 



d , , 




d 


dt^^^'l 


N 


~^ It 



(21) 

under the integral in eq. p7)) . The quantity E^^^ with this replacement will be 
called Ef,. The over-bared quantities refer to the bosonic scattering data. The re- 
placement eq. (|21[) must, of course, be accompanied by the boson Feynman dia- 
gram i?B so that the total vacuum polarization energy becomes 

-Evac = i?5 + £^fd' + J (22) 

where the superscript indicates the inclusion of the counterterm contributions. Each 
of the three terms on the right hand side of equation (|22p is ultra-violet finite by 
itself. The advantage of eq. ([22l) and the replacement eq. ([2T|) is now obvious: Instead 
of fermionic contributions up to order A'^ = 4, we only need to compute second 
order fermionic and bosonic Feynman diagrams and terms in the corresponding 
Born series. 
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4. Numerical Results for the Vacuum Polarization Energy 

Numerical results as well as dimensionfuU parameters are measured in appropriate 
units of the (perturbative) fermion mass m. 

We focus on the contribution of fermion fluctuations to the vacuum polarization 
energy because it dominates the boson counterpart by a factor proportional to the 
number of internal degrees of freedom, e.g. Nq, the number of colors. In this scenario 
we are now prepared to compute the vacuum polarization energy of a prescribed 
string configuration. 



4.1. Variational Ansdtze 

Despite of the simplification in eq. ([22l) , the numerical computation is still expensive. 
The scattering data are extracted from a multi-channel problem and for the final 
result to be reliable a huge number of partial wave must be included. This numerical 
effort restricts the number of variational parameters that can be used to characterize 
the profile functions. We have already introduced the strength parameter ^i. In 
addition, we introduce three scale parameters vjw and via the ansatze 



/ff(p) = l-e-^, fa{p) = l-e UgJ , ^(p) = 6 



1 



(23) 



The scale parameterizes the shape of the gauge profile. As explained above, this 
shape and thus should not be observable. The other specifics of the profiles are 
chosen to keep E^i regular. 

We have also considered an exponential parameterization for the gauge field 



/g(p) = 1- 1 + ^ exp , (24) 

V wgJ \ wgJ 

which yields a slightly better agreement with the original Nielsen-Olesen profiles 
that minimize -Ed for = 7r/2. No significant difference in -Evac was found between 
these ansatze. 



4.2. Gauge Invariance 

We check gauge invariance by varying the shape of the gauge profile, (,{p). A typical 
result is shown in table [TJ As expected, the individual contributions to i?vac depend 

Table 1. Numerical results for the various con- 
tributions I I22I I to the fermion vacuum polariza- 
tion energy in the minimal subtraction scheme. 





^FD 


Es 


^rcn. 




0.5 


-0.2515 


0.3489 


0.0046 


0.1020 


1.0 


-0.0655 


0.1606 


0.0032 


0.0983 


2.0 


-0.0358 


0.1294 


0.0038 


0.0974 


3.0 


-0.0320 


0.1235 


0.0056 


0.0971 


4.0 


-0.0302 


0.1193 


0.0080 


0.0971 
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strongly on w^. However, these changes essentially compensate each other. Numer- 
ically the most cumbersome part of the calculation is Eg . From various numerical 
considerations (change of extrapolation scheme for partial wave sum, modification 
of momentum integration grid, etc.) its numerical accuracy is estimated to be at 
the 1% level. Within that range -Evac is independent of w^, thus verifying gauge 
invariance. 

The above results are obtained in the MS renormalization scheme, which essen- 
tially omits the non-divergent parts of the Feynman diagrams. Any other scheme 
merely differs by manifestly gauge invariant (finite) counterterms. 

4.3. On-Shell Renormalization 

With the above mentioned choice of units, the dependence of E^ac on the model 
parameters factorizes in the MS scheme which simplifies the computation because 
this dependence can easily be traced from E^. However, for physically meaningful 
results we need to impose renormalization conditions that correspond to a particle 
interpretation, inducing a mild parameter dependence in Eqt- To be specific we 
consider the so-called on-shell scheme in which the coefficients of the four allowed 
counterterms are determined such that 

• the tadpole graph vanishes 

• the Higgs mass remains unchanged 

• the normalization of Higgs particle remains unchanged 

• and the normalization of vector meson remains unchanged 

in the presence of fermionic quantum corrections. Note that the vector meson 
mass Mw is not fixed by these conditions and thus will be a prediction that includes 
quantum corrections. Hence we tune the gauge coupling to reproduce the physical 
value for Mw- Typical results for the vacuum polarization energy per unit length of 
the string are shown in figure [U as functions of the variational parameters. Except 
for narrow string configurations dominated by the Higgs field, the vacuum polar- 
ization energy turns out to be positive. Therefore, fermionic quantum corrections 
do not provide any sensible binding and no stable uncharged string is found for the 
physically motivated parameters, eq. ([7]), for which E^i dominates the total energy. 
Yet, Ec\ deceases quickly with increasing Yukawa coupling / and some stability is 
indeed seen for large / and narrow strings. Unfortunately, in this regime the restric- 
tion to one fermion loop in the vacuum polarization energy is unreliable because of 
the occurrence of the Landau ghostP^. 

5. Charged Strings 

Cosmic strings induce many fermionic bound state levels (whose energies are de- 
noted by tj) for the two-dimensional scattering problem. For = 7r/2 there even 
exists an exact zero mod4^. In the three dimensional problem these bound states 
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acquire a longitudinal momentum for the motion along the symmetry axis and their 
energies become 



-Pn 



with 



Pn = -r 



nir 

T 



(25) 



Here L is the length of the string. In leading order of the L — )■ cx) limit the 
sum over the discrete longitudinal momentum turns into a continuum integral, 
— > J dp. To minimize the bound state contribution a chemical potential /i 
with max(|ej|) < < m is introduced and all levels with Ei[p) < /i are populated. 
This procedure defines a Fermi momentum for each level, pf (/i) = y/ ii^ ~ ef which 
enters the total charge per unit length of the string 

Y-pf(A^) 



(26) 



This relation can be inverted to give /i — fJ.{Q) and thus pf = pf (Q). From this 
the binding energy (per unit length) for a prescribed charge 

i^bind(Q) - - V / dp 



p^ 



(27) 



is computed relative to an equal number of free fermions that have energy m each. 
Figure [5] shows the fermion contribution to the binding energy, -Evac + £'bind('3)- For 
a given configuration the graph terminates at the point when all available bound 
state levels are occupied. For small charges narrow strings are favorable while the 
binding energy of strings with larger widths decreases more quickly as Q increases. 
Surprisingly, the envelope along which -Evac + E'bind is minimal forms a straight line. 
Extrapolating this line to Q = indicates that the fermion vacuum polarization 
energy should (approximately) vanish. This extrapolation circumvents the Landau 
ghost problem. 



"T 1 ' 

: 0.5 = 0.5 

0.5 = 2.0 

2.0 = 0.5 

2.0 = 2.0 



0,2 0,3 
J; /It 



Fig. 1. Fermion vacuum polarization energy in the on-shell renormalization scheme. 
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To finally decide on dynamical stability, the classical energy must be included. 
To this end we scan through several hundred configurations characterized by the 
variational ansatze, eq. (^5)). We label them by s = 1, 2, . . . and compute their total 
binding energy 



is) 



Nr. 



for a given charge. If 



EtotiQ) :-min, E^'JiQ) 



< 



(28) 



(29) 



a stable configuration is constructed. Figure [3]shows Etot as a function of charge for 
various values of the Yukawa coupling constant, i.e. the mass of a non-interacting 
fermion. For / ss 1.6 the classical and fermion energies essentially cancel each other 
and leave -Etot roughly charge independen10. Bound objects are observed by fur- 
ther increasing the Yukawa coupling to about / « 1.7, which corresponds to a 
heavy fermion mass which is less than twice the top quark mass. We find that the 
minimizing configurations have ~ 0, i.e. they are dominated by the Higgs field. 



6. Conclusion 

We have presented and discussed the formalism to compute the fermion contribution 
to the vacuum polarization energy per unit length of an infinitely long straight 
string in a simplified version of the electroweak standard model. Our approach is 



'^The exhibited dependence at small Q is artificial because very narrow strings have not been 
considered to avoid the Landau ghost inconsistency. 




Fig. 2. (Color online) Fermion contribution to the energy of a charged string. 
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based on the interface formalism, for which the analytical properties of scattering 
data are essential. We have also seen that a particular subset of gauge choices 
circumvents obstacles that in a naive treatment arise from the non-trivial structure 
of the string configuration at spatial infinity. Numerically we have found that the 
vacuum polarization is small and positive in the regime in which the one-fermion 
loop approximation is reliable. Hence, there is no quantum stabilization of the string. 
However, we have seen that a heavy fermion doublet can stabilize a nontrivial string 
background for a non-zero fixed charge per unit length. The resulting configuration 
is dominated by the Higgs field. Since any additional variational degree of freedom 
can only lower the total energy, the embedding of this configuration in the full 
standard model, with the U{1) gauge field included, will also yield a bound object 
as long as mixing between this heavy and the standard model fermions can be 
ignored. We see binding set in at m sa 300 GeV, which is still within the range 
of energy scales at which the standard model is expected to provide an effective 
description of the relevant physics, and also within the range to be probed at the 
LHC. Light fermions would contribute only weakly to the binding of the string, since 
their Yukawa couplings are small. As a result, we can add them to our model, e.g. 
to accommodate anomaly cancellation, without significantly changing the result. 
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